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Application of a sufficiently strong par- 
allel magnetic field B» > B c produces a 
soliton-lattice (SL) ground state in a double- 
layer quantum Hall system. We calculate the 
ground-state properties of the SL state as a 
function of Bu for total filling factor z/p = 1, 
and obtain the total energy, anisotropic SL 
stiffness, Kosterlitz-Thouless melting tem- 
perature, and SL magnetization. The SL 
magnetization might be experimentally mea- 
surable, and the magnetic susceptibility di- 
verges as \Bh — -Bel -1 . 

Keywords: quantum Hall, solitons, magne- 
tization, Kosterlitz-Thouless 

I. INTRODUCTION 

At sufficiently small layer separation, a 
double-layer quantum Hall (2LQH) system 
is an unusual quantum itinerant fer- 
romagnet. |||4| The 2LQH system can be 
mapped to an equivalent spin- 1/2 system by 
equating "up" ("down") pseudospins with 
electrons in the upper (lower) layer. || For 
any finite layer separation d, the itinerant fer- 
romagnet has an XY symmetry, so that the 
orientation of a pseudospin at location r is 
specified by its angle 0(r) in the xy plane. 

Murphy et al. have investigated the effect 
of an in-plane magnetic field B\\ on 2LQH sys- 
tems and find evidence for a phase transition 
between two competing QH ground states at 
a critical value B\\ = B c . || These two ground 
states have been explained theoretically 
by showing that B\\ produces a rotating Zee- 
man field seen by the pseudospins. This gives 



a pseudospin contribution to the ground- 
state energy of the Pokrovsky-Talapov (PT) 
form, §j7[ 
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where Q = 2ttz x B|id/^>0 defines the paral- 
lel magnetic-field wave vector, 9(r) = 9(r) + 
Q • r, t = t e~ Q e ^\/A&\Vi is the tunnel- 
ing energy (t is the tunneling energy when 
Q = 0) B an d Ps = Pe{^v\i>2) is the pseu- 
dospin stiffness [pE is the interlayer exchange 
stiffness when V\ — 1/2 — 1/2) in the Hartree- 
Fock Approximation (HFA), Uj is the fill- 
ing factor of layer j, and the energy is mea- 
sured relative to the ground-state energy for 
B\\ = 0. Note that by adjusting the front 
and back gate voltages of the sample, V\ and 
vi may be varied (with z/p = v\ + V2 = 1), 
thereby allowing t and p s to be adjusted. 

For small B\\, Eq. ( |1 . 1| ) is minimized by 
9(r) = 0. This is the commensurate (C) 
ground state. For all finite Bu > B c , the 
pseudospin polarization has broken transla- 
tional symmetry, and a soliton-lattice (SL) 
state results. For large Bu the pseudospins 
behave (almost) as if t = 0. This work fo- 
cusses on calculating the ground-state prop- 
erties of the SL state, for all Bu > B c . Inter- 
estingly, it is not necessary to solve for the 
form of 6(r) in order to calculate the total 
energy of the system. |6|j7],f| 

Minimizing E with respect to 9 gives the 
2D sine-Gordon equ ation (S GE), £ 2 V 2 # = 

sin#, where £/£ = J2np s /t. We shall give 



numerical values of our results for a hypo- 
thetical GaAs 2LQH sample with total den- 
sity 1.1 x 10 n cm -2 , layer separation d = 21.1 
nm, and tunneling energy /j = 0.1 meV. Such 
a sample would have I ~ 11.8 nm, d/£ = 1.8, 
and hu c ~ 8 meV for z/p = 1, and pe ~ 0.08 
meV and £ ~ 26.5 nm in the HFA. 

II. SOLITON-LATTICE STATE 

We investigate solutions of the SGE of the 
form 9(r) — B [e^ • (r — r )], where §x is some 
unit vector in the xy plane, so that £, 2 df6 = 
sin 9, which is closely analogous to the equa- 
tion of motion of a pendulum. The conserved 
quantity analogous to the total energy of a 
pendulum is 2c 2 = (l/2)£ 2 (<9i#) 2 -(l-cos#). 
Defining (3 = 6/2 leads to the equation 



£<9i/3 = ±Vc 2 + sin 2 /?. 



(2-1) 



When c = 0, Eq. (Q gives 9 ss (r) = 
4arctanexp[ei-(r — ro)/£], which represents a 
single soliton of width £, corresponding to the 
motion of a pendulum that completes a sin- 
gle revolution in infinite time. The energy per 
unit length of a single soliton relative to the C 
ground-state energy may be computed from 
Eq. (0) as E ss /L 2 = p s (8/f - 2^ • Q), 
where L 2 is the sample length perpendicu- 
lar to ei. The lowest energy solitons have 
B\ = Q and their energy vanishes when 
Q = Q c = 4/(7r£). Thus for Q > Q c it is en- 
ergetically favorable to create solitons. The 
solitons are weakly repulsive, and form a SL. 
Analogous soliton effects also occur in long 
Josephson junctions. 
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The SL spacing L s may be determined 
by noting that over one period of the SL, (3 
changes by it. Thus 
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where we have used Eq. fl2.1|) , defined 77 = 
l/Vc 2 + 1, and where X(^) is the complete 



elliptic integral of the first kind. ]TT[ We de- 
fine the SL wave vector Q s = (27r/L s )ei, so 
that Eq. (|2.2|) becomes 



(2.3) 



Q s /Q c = (n/2) 2 /{ V K(r))}. 



Note that r\ — > corresponds to Q s — >■ Q — ► 
00, whereas 77 — > 1 corresponds to the C-SL 
transition, where Q ^ Q c and Q s ■— >• 0. 

The energy per unit area is obtained by 
expressing Eq. (|1 . 1|) as an integral over (3 [cf. 
Eq. (|2.2| )1, which gives 
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where E{rj) is the complete elliptic integral 
of the second kind |l l] and L\L 2 is the sam- 
ple area. The value of Q s that minimizes 
the energy per unit area is found by dif- 
ferentiating Eq. (2.4) with respect to Q s , 
holding Q constant. Using the identity JTTJ 
dE/dr] = [E(i]) — K(j])]/r), one obtains 



Q/Q c = E(fj)/fj. 



(2.5) 



Equations ( |2.3| ) and (|2.5|) together determine 
the equilibrium SL wave vector Q S (Q) that 
minimizes the energy. For Q/Q c — > 00, 
Qs/Q « 1 - (l/2)[(vr/4)(Q c /g)] 4 . For Q 
Qc Qs/Qc ~ (^V 2 )/ 111 ! 1 / 6 ) asymptotically, 
HJl2| where e = Q/Q c — 1. If it were pos- 
sible to achieve e ~ 10~ 2 (e.g., by gating 
the sample to tune Q c ), surface acoustic wave 
(SAW) techniques might detect the SL when 
L s matched the SAW wavelength. The gen- 
eral SL ground-state solution for 9 is given 
by sin^sL - 7r)/2] = sn[ei ■ (r - r )/r/£,r/], 
where sn is the sine-amplitude Jacobian func- 
tion. [|TIJ 



III. STIFFNESSES 

The elastic constants of the soliton lattice 
are given by the stiffness tensor 
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which we calculate here for all Q > Q c using 
Eqs. (2.4) and (|2.5|). Our results agree with 
those obtained in Ref. |12| for Q — > Q c . 

Using the results of Sec. |J and the iden- 
tity dK/drj = E( V )/[rj(l-rj 2 )}-K(rj)/rj, 
the compressional elastic constant A n is 
found to be 



89 



(3.4) 



dQ 
dQ s 



1 3 / 7T Qc 

2 \4 Q 



(2/7r 2 )eln 2 



l/e), Q/g c 



(3.2 



As expected, A n — > p, 
oo; K n -> for Q/Q_ 
short-range repulsions between the solitons. 
The shear elastic constant A 22 is 



in the limit Q/Q c — > 
— > 1 because of the 



where we have integrated by parts and made 
use of the SGE. The energy SE is minimized 
by choosing 89 from among the eigenvalues 
of bracketed Schrodinger-like operator in Eq. 
(p.4|). For Q = x, 6*o = ^(x) and we may 
write S9(r) oc exp(iq y y)89(x). When 9 (r) = 
9sl{ y ), the eigenvalue equation for brack- 
ooeted equation in Eq. (|3.4|) becomes Lame's 
1 equation after a simple rescaling of x. JTH 
The relevant "vortex oscillation" solutions to 
Lame's equation follow from Ref. [10[] and 
Eqs. ( |2.3| ) and Q2.5[ ) in the long- wavelength 
limit: 8EjL\L 2 
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Like An, approaches A 22 — > p s for Q/Q c — > 
oo. However, A 2 2 diverges for Q —> Q c . 

The SL phase of the PT model can un- 
dergo a KT transition, |7| with a transition 
temperature /csTkt ~ (^/^VKuK^. The 
KT transition is probably most easily mea- 
sured in devices with oppositely directed cur- 
rents in each layer. M However, because the 
SL dislocations are electrically charged, the 
KT transition might increase the longitudinal 
resistivity even in devices without separately 
contacted layers, due to the proliferation of 
unbound charged dislocations above T K t- 



We now calculate the stiffness tensor A, 



relevant to the sound velocities of the SL by 
writing 9(r) = 9q(t) + 89(r), where 9q is the 
ground-state solution of the SGE which min- 
imizes the PT energy ( |1 . 1|) and 59 is a small 
deviation of 9 from 9 . To quadratic order in 
9, the change in the PT energy ( |1 . 1| ) of the 
2LQH system is 



(1/2) (A n g 2 + p s g 2 ), where 
An is equal to the compressional stiffness 
(p|). Thus An = A n , while A 22 = p s is 
independent of Q. 

IV. MAGNETIZATION AND 
MAGNETIC SUSCEPTIBILITY 

The SL makes a small but possibly mea- 
surable contribution to the magnetization of 
the 2LQH system. The parallel-field magne- 
tization per unit area is 
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where the equilibrium value of E/L\L 2 , 
found by using Eqs. ( |2.3|) and ( j2.5|) in Eq. 
(2.4), and we have used the results of Sec. [TI|. 
M = 2Ttp s Q c d/(j>o sets the scale of the mag- 
netization density, which we plot in Fig. 1. 
It is useful to compare the total SL magne- 
tization M LiL 2 to the scale of the Landau 
diamagnetism in a u = 1 QH system: 
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where p^ = eh/2m*c is the effective Bohr 
magneton in GaAs, and m* « 0.067m e is the 
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effective mass. This shows that in the HFA, 
the SL magnetization is roughly an order of 
magnitude smaller than the Landau diamag- 
netism, which has been measured in GaAs 



single-layer heterostructures. ]T3|] The torque 
on a tilted 2LQH sample has both a Landau- 
diamagnetic component t± = M±B\\ and a 
SL component tii = M\\B±. The smallness 
Mm is somewhat compensated by the size of 
B± in the expression for rii. 

The parallel-field magnetic susceptibility 
X\\ = dM/dBu is given by 
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Xo(dQJdQ - 1) = Xo(p s /K n - 1) (4.3) 
(3/2)[(tt/4)(Q c /Q)] 4 , Q/Q c ^oc 
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where xo = (27rd/0o) As ~ 5 x 10 m. Near 
the C-SL transition, x\\ ~ V e > with logarith- 
mic corrections. It might not prove practical 
to make AC measurements of x\\ because of 
substantial eddy-currents in the QH regime. 
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FIG. 1. SL contribution to the magnetization 
density (dots), which drops precipitously for 
Q > Q c - The C-phase result (solid line) holds 
for Q < Q c . 
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